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Abstract 

We investigate the thermodynamics of an integrable spin ladder model which 
possesses a free parameter besides the rung and leg couplings. This model is exactly 
solvable by means of the Bethe Ansatz and exhibits a phase transition between a 
gapped and a gapless spin excitation spectrum. The magnetic susceptibility is 
obtained numerically and its dependence on the anisotropy parameter is determined. 
A comparison between the spin gap obtained from the susceptibility curve and the 
one obtained from the Bethe Ansatz Equations is made and a good agreement found. 
A connection with the compounds KCuCls dill, CM2(C5-ffi2^2)2C'^4 [BlllllSj and 
{C^Hi2N)2CuBr4^ [B] in the strong coupling regime is made and our results for the 
magnetic susceptibility fit the experimental data remarkably well. 

PACS: 75.10. Jm, 71.10.Fd, 03.65.Fd 



1 



1 Introduction 



The study of spin ladders, both from an experimental as well as from a theoretical point 
of view has gained much attention in the last few years and has by now a literature of 
its own. On the one hand with the discovery of high-Tc superconductivity in doped 
cuprate materials |7j, a tremendous effort has been made to understand the physics of 
these compounds, in which the Cu-0 ions build up a two-dimensional lattice. This 
task is made more difficult by the lack of exactly integrable models as compared to 1-D 
where an abundance of such models has been known for decades. On the other hand 
spin ladders are prototypes of quasi one-dimensional systems and as such represent an 
excellent proving ground for studying the transition from the more amenable and better 
understood physics of one-dimensional quantum systems to two-dimensional behavior. 
They exhibit antiferromagnetic (AF) ground states, where at each site of the lattice (in 
the undoped state) an electron occupies one of two spin states As such they can be 
reasonably well approximated by the archetypical AF Heisenberg model or some suitable 
generalization thereof. In one dimension the Heisenberg model is exactly solvable via the 
Bethe Ansatz and from the solution it is known that in the AF regime, up to a given 
threshold of the spin exchange interaction, the elementary spin excitations are gapless 
(beyond the threshold one has a Neel-like groundstate with gapped excitations) PI ITU]. 
The existence of a spin gap is critical for the occurence of superconductivity under doping, 
whereby the holes introduced through doping undergo a Bose-Einstein condensation. By 
introducing the concept of a ladder model this apparent contradiction is resolved, since 
these systems allow for the formation of singlet states along the rungs which in turn 
are responsible for the appearance of a spin gap. (Strictly speaking the singlet states 
can only form when there is an even number of legs fUEl)- Recently, with the rapid 
progress being made in nano-engineering, many different species of ladder compounds 
such as SrCu203, Lai_xSrxCu02.5, Srn_xC axCu2i04^i and other organic compounds j2] 
have been investigated in this context. This theory is now supported by a substantial 
body of experimental evidence. 

However, in contradistinction to its one-dimensional counterpart, the two-dimensional 
Heisenberg ladder model cannot be solved exactly. So in the search for exactly integrable 
models which might give one a better insight into the physics of ladder systems, many 
authors have considered generalized models which incorporate additional interaction terms 
without violating integrability [T^14pi5pi6pi7pi8pi9p20p2ip22p23|l^ . Notwithstanding the 
fact that these models exhibit realistic physical properties such as the existence of a spin 
gap PHO^ and the magnetization plateaus at fractional values of the total magnetization 

, up to our knowledge none has been used to predict physical quantities that could be 
compared directly with experimental data. 

The main goal of the present work is to show that integrable ladder models can be used 
to model experimental data. In particular, we discuss a generalized integrable spin ladder 
with incorporates an additional anisotropy parameter besides the rung and leg interac- 
tions and show that the magnetic susceptibility as obtained from this model as a function 
of temperature fits well the experimental data for the ladder compounds KCuCl^ DDIS, 
Cu2{Cr,Hi2N2)2Ch mmSl and {C^Hi2N)2CuBri 6j. The model is exactly solvable by 
the Bethe Ansatz and reduces to the model introduced by Wang for some special 
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limit of this extra parameter. This additional free parameter, which we call t, brings 
about an anisotropy-dependent critical line through an explicit functional dependence of 
the excitation gap. The thermodynamics of the model is investigated and the magnetic 
susceptibility curve as a function of the temperature is obtained. The influence of this 
anisotropy parameter on the thermodynamics of the model is also determined. A compar- 
ison between the spin gap obtained from the susceptibility curve and that one obtained 
from the Bethe ansatz equations is performed and a good agreement is found. 

The paper is divided as follows: In section 2 we describe the model and its Bethe 
Ansatz solution. Section 3 discusses the thermodynamical properties of the model and 
is also devoted to a comparison between the experimental results and our findings. We 
conclude the paper with a brief summary. 



2 The model 

Let us begin by presenting the integrable anisotropic spin ladder model in the presence 
of an external magnetic field h |23], whose Hamiltonian reads 

where 

Periodic boundary conditions are imposed along the legs. As usual cfj and Tj are Pauli 
matrices acting on site j of the upper and lower legs respectively, Jr{Ji) is the strength of 
the rung coupling (leg coupling) and t is a free parameter which introduces an anisotropy 
in the leg and in the interchain interactions. Throughout this paper L stands for the 
number of rungs (the length of the ladder). By setting t — 1 in Eq. one recovers 
the isotropic model of Wang [23] based on the S'[/(4) symmetry (strictly speaking Wang's 
model is S'f/(4)-invariant only in the absence of the rung interactions). The Hamiltonian 
is invariant under the interchange of the legs, i.e. (?j ^ fj. Moreover, under spin inversion 
on both leg spaces it is also invariant under the exchange t <-> t~^. The energy eigenvalues 
of ((H) are given by 

|-!:(^-^^)-0-^f)^-^w--3) 

where the A/s are solutions of the three-level nested Bethe ansatz equations 

\\ + VSy* A, -M + i 11 A, - - z/2 
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For the sake of convenience we take from here on J; = 1. From an analysis of these 
equations it follows that for > 1 + ^(t + j) the reference state becomes the ground 
state and any excitation is gapped. In this region the ground state is given by a product 
of rung singlets, indicating a dimerization along the rungs. The energy gap A can be 
calculated using the exact Bethe ansatz solution and has the form 



A = 2 X - 1 



h 1, 1. 

t + - 

2 2^ t' 



(4) 



By solving A = for Jr we find the critical value J^^ = 1 + + j), indicating the critical 
line at which the quantum phase transition from the dimerized phase to the gapless phase 
occurs. This result is depicted in Fig. 1. 

Notice that by choosing t real one has Jr"^ > 2 which corresponds to the strong coupling 
limit. We therefore expect our anisotropic model to be a good candidate for describing 
some organic ladder compounds, such as KCuCl^, Cu2{C^Hi2N2)2Cli, TlCuCl^ and 
{C^Hi2N)2CuBri. 















gapped 










Wang's point 




gapless 







Figure 1: The phase diagram of (0) in the (J^. x t)-plane in units of Ji. The arrow 
indicates Wang's point (t=l). The curve Jr^ = 1 + {t + l/t)/2 divides the gaped from the 
gapless phase. 



3 Comparison with experimental data 

Our main goal is to calculate the magnetic susceptibility x of our exactly integrable ladder 
Hamiltonian and compare it to experimental data on the compounds mentioned before. 
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In order to do this we block-diagonalize the Hamiltonian Eq. for a sequence of finite 
chains (up to 6 rungs or 12 sites) and for different values of the couphng and anisotropy 
t. With these spectra we generate a sequence of finite-size susceptibihties x^^^ as a function 
of the temperature T and then apply a standard Bulirsch-Stoer acceleration technique 
P7] to the sequence in order to determine the susceptibility in the thermodynamic limit 
L — >■ oo. An example for = 7 and t = 1 is presented in Fig. 2. In order to compare 



J = 7 and t = 1 




Figure 2: The finite-size susceptibility per site for the case Jr = 7 and t = 1 for varying 
number of rungs L and the extrapolated curve. 



our results with experimental data available on a series of compounds, a few remarks on 
the behavior of the susceptibility as a function of Jr and t are necessary, as well as the 
procedure used for choosing the values of our parameters. First we observe that for a fixed 
value of t there is a 'smoothing out' of the susceptibility curve for increasing Jj. while for 
fixed Jr the magnetic susceptibility increases with decreasing t, as can be seen in Fig. 3a. 
The second point is that in all cases the susceptibility presents an exponential decay for 
low temperature T << A* where A* is the spin gap of the system 

A* 

e T 



T 

which is in agreement with the result of Troyer et. al. j2H] for the Heisenberg spin ladder 
model. By linearizing Eq. ^ a numerical value for the spin gap (A*) can be found. 
When compared to the exact expression of the spin gap A obtained from the Bethe 
Ansatz Equations (j3]) at T = an excellent agreement is found as can be seen in Table 
1. The linearized curves are depicted in Fig. 3b. 

We remark that the behavior of the magnetic susceptibility for the critical value of 
the couplings (J^^^) is of the form 

X- ~ (6) 



which indicates a typical quantum critical behavior. This had already been predicted 
in PHOH] for the isotropic case as illustrated in Fig. 4 for different values of t. By 
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(a) 



(b) 




Figure 3: a) The magnetic susceptibility (x) as a function of the temperature (T) for 
different values of the coupling (Jr) and anisotropy (t). (b) A logarithmic plot of the 
susceptibility (x) as a function of the inverse of the temperature (1/T), from which the 
spin gap (A*) can be obtained. 
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0.30 


10.37 


10.36 




0.20 


8.80 


8.81 




1.00 


20.00 


19.88 


12 


0.30 


18.37 


18.32 




0.20 


16.80 


16.76 



Table 1: Spin gap A* obtained from the linearization of Eq. (5) and Fig. 36 compared to 
the analytical result A obtained from the Bethe Ansatz Equations (4) for different values 
of Jr and t. 




t=1 .00 , J '=2.00 

t=0.75 , J,'=2.04 

t=0.50 , J,'=2.25 



■ 0.0 2.5 5.0 

T 

Figure 4: Magnetic susceptibility (x) versus temperature (T) for different values of the 
coupling Jr in units of Ji. These graphs indicate a typical quantum critical behavior. 

setting our parameters accordingly we were able to obtain the spin gap for some ladder 
compounds such as KCuCk [Uig, Cu2{C5H^2N2)2Ch HHE], {C^Hi2N)2CuBr^ which 
were previously studied by different authors using different methods. A comparison of 
our choice of parameters with the experimental data A*^^^ [211 for the above mentioned 
compounds is depicted in Table 2. For the sake of clarity we present the results on three 
different compounds separately and refer the reader to the quoted works for more details 
on the experimental setup. 

3.1 KCuCh 

The magnetic properties of this system were investigated in P3I2] where a non-magnetic 
ground state and a spin gap in the excitation spectrum were reported. The crystal struc- 
ture of this compound is moniclinic and the double chain feature of the system arises from 
two edge-sharing chains of CuCIq octahedra (see [Hj for more details on the geometry of 
this compound). 
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compound 


Jr 


t 


A 


A* 




KCuCh 


4.0 


0.3254 


32.00 


36.56 


31.10 


Cu2iC,Hl2N2)2Ch 


5.5 


0.2344 


10.80 


13.41 


10.80 


{C5Hi2N)2CuBu 


4.0 


0.3100 


9.85 


9.90 


9.50 



Table 2: The values of parameters for different compounds. The experimental data were 
obtained from Ref. 



The experimental data 1| and exact diagonalization results are depicted in Fig. 5. 



J= 12.3 J = 49.2 and 1 = 0.3575 




T 



Figure 5: Magnetic susceptibility versus temperature. The full line is the result from 
extrapolation of finite-size data obtained from exact diagonalization of chains up to 12 
spins. Circles are experimental data obtained by Nakamura and Okamoto pp. The inset 
shows the low-temperature regime. 

The spin gap we used was estimated experimentally to be ~ 31. 10/^ [2^. The full 
line was obtained using the following exchange couplings: J/ = 12. 3i^', Jr = 49.2K and 
t = 0.3575. We observe that the low-temperature regime shows a better agreement 
between theory and experiment. The deviation in the high-temperature regime might 
arise from the existence, in the compound, of a ferromagnetic coupling which takes place 
between the 3d orbitals of Cu^"^ ions mediated by the CI anions [30^, a coupling which is 
not present in our model. 

3.2 Cu2{C^HuN2)2Ch 

This compound have been studied by Chaboussant et al The ladder structure 

arises from the stacking of Cu-Cu binuclear units (each Cu is to be seen as belonging 
to one leg). For magnetic fields below h = 7.5T the one-dimensional ground state is a 
valence-bond singlet structure on each rung and the system is gapped. Above this lower 
critical field, the first excited state (a triplet) becomes the ground state and the system has 
a non-zero magnetization. There are basically three different exchange paths between Cu 
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ions: the intramolecular exchange constant Jr between Cu-Cu ions takes place through 
two CI ions while the intermolecular Ji involves in-plane unpaired electron densities of 
Cu ions of different molecules. The third intermolecular exchange introduces a frustration 
via next-nearest neighbor interaction and occurs through hydrogen bonds between the 
CI and N ions. This interaction is however significantly smaller than J^. and is taken 
altogether equal to zero in our model. The low-temperature gap has been determined 
through susceptibility measurements and NMR techniques to be A = 10.8 ± 0.6K with 
Ji = 2.4K and Jr = 13. 2K The same gap can be obtained in our model with the choices 
Ji = 2AK, Jr = 13.2K and t = 0.2344. The agreement between the low-temperature 
experimental data and our results is excellent, as can be seen in Fig. 6, and remains good 
for the whole temperature range. 

J|^2.4 J= 13.2 and 0.2344 

1,0 
0,8 
0,6 
0,4 
0,2 



"'"0 10 20 30 40 

T 

Figure 6: Normalized magnetic susceptibility versus temperature. The full line is the 
result of exact diagonalization followed by extrapolation of finite-size data. Again we used 
chains of up to 12 spins. The circles are the experimental data obtained by Chaboussant 
et al j4j. The inset depicts the low-temperature regime. 




3.3 {C5HuN)2CuBr4 

As in the previous compound, here we also have a binuclear Cu-Cu stacked structure. 
The intramolecular exchange coupling is mediated through an orbital overlap of Br ions 
adjacent to the Cu sites. The intermolecular exchange is also mediated by Br ions a bit 
further apart (8.597A in comparison to intramolecular distances of 6.934A). A diagonal 
frustration exchange should be present but it is expected to be weak 6J. The exchange 
couplings used were Ji = 4.0i^, Jr = 16. Oi^ and t = 0.31i^r. Again the fit is better for low 
temperatures, with an energy gap of 9.8 K (9.5 K experimentally [U]) and an exponential 
drop in the magnetic susceptibility. The data is depicted in Fig. 7. 
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Figure 7: Magnetic susceptibility versus temperature. Circles represent experimental 
data obtained by Watson et. al. The solid line is the result obtained via exact diago- 
nalization as explained in the previous captions. The inset depicts the low-temperature 
regime. 

4 Summary 

We presented a numerical analysis of thermodynamic properties of an integrable anisotropic 
spin ladder model which has, besides the rung and leg couplings, an extra free anisotropy 
parameter. From the Bethe Ansatz equations as well as the exact expression for the 
energy spectrum we showed that the model exhibits an anisotropy- dependent spin 
gap. The critical line and phase diagram were obtained. The magnetic susceptibility was 
investigated and good agreement was found between the gap obtained from the Bethe 
ansatz equations and the gap obtained from the magnetic susceptibility. A connection 
with some strong coupling compounds was discussed and our results reproduce well ex- 
perimental data for these compounds in strong coupling limit. The anisotropy parameter 
in our model has no obvious physical interpretation. However, cyclic four-spin exchange 
terms are known to affect the excitation spectrum of spin ladders and the importance 
of this type of interaction for appropriately describing the properties of strong-coupling 
ladders has also been established [H2]. As shown in this work the anisotropy parameter 
allows us to tune our model to different physical systems of interest, thus presenting a 
unified scenario for strong-coupling ladders. 

This work has been financially supported by the brazilian federal agency CNPq (APT, 
AF and AIM) and the state agency FAPESP (AIM). 
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